Exchange energy dominated by large orbital spin-currents in 6-Pu 
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The electronic structure ol the anomalous (5-phase of Pu is analyzed by a general and exact 
reformulation of the exchange energy of the /-shell. It is found that the dominating contribution 
to the exchange energy is a polarization of orbital spin-currents that preserves the time reversal 
symmetry, hence a non-magnetic solution in accordance with experiments. The analysis brings a 
unifying picture of the role of exchange in the 5/-shell with its relatively strong spin-orbit coupling. 



^— ■») ' The results are in good accordance with recent measurements of the branching ratio for the d to f 

^— ■>) ' transition in the actinides. 

(N '• 

The electronic structure of the actinides stands out among the elements of the periodic table as most intricate. Their 
, ■ 5/ states form narrow bands with a comparatively large spin-orbit coupling. Rather subtle changes in the electronic 
, structure leads to different ground states, and in various actinide compounds the 5/ show itinerant, localised as 
well as heavy-fermion behaviour and are responsible for both magnetic order and superconductivity. In recent years 
■ this complexity has been exemplified by various attempts to theoretically understand the phase diagram of Pu, and 
especially the formation of its high temperature, large volume, highly anomalous (5-phase [H, [3jSflM00, S|- As a 
first progress, it was observed that the stability of this phase can be understood if allowed for spin polarisation (SP) 

' i>l2|- However, the existence of any magnetic moment are in contrast to a large amount of experimental observations 
I , 9|. Most recently it was pointed out that orbital polarisation (OP) plays a major role in stabilising this magnetic 

' solution 0. Then, by utilizing the so-called LDA+U approach, i.e. a local density approximation to density functional 

1 t ■ theory with an added extra local Hartree-Fock (HF) term. When using a HF term in its most general form allowing 

' for off-diagonal spin coupling [1, Q , it was observed that a non-magnetic large volume phase could be stabilised in 
^ , a somewhat counter intuitive way - by increasing the exchange energy the moments vanish. In the analysis of these 
^ ' calculations, it was pointed out that the solution was closer to a jj coupling scheme than the usual LS coupling and 



that the configuration was closer to a than the expected /^. An observation of such a co nfig uration is again 
opposing many experimental findings, although some recent experiments verify a jj-like coupling l3, lH. Meanwhile 



o . 

O . there have also been various dynamical mean field theoretical (DMFT) calculations, which also lead to a high volume 
non-magnetic state, at least when allowing for the off-diagonal spin coupling, which they attribute to a localization 

, of the 5/ states. There are also a model of mixed valency, where part of the 5/ states are localised while the rest 
J> ' are itinerant With four localised and one itinerant a large volume phase is calculated to be stabilized with an 

, electronic structure close to what is observed by photo-emission experiments [l2j . 

' The purpose of the present Letter is to explore the relationship of the SP-I-OP approach with that of the general 
, LDA-I-U approach with its more flexible exchange interaction, and to explain how the non-magnetic state is stabilised 
by the latter. This is accomplished be re-expressing the general exchange interaction of the LDA-I-U approach as a 
sum of interacting multi-poles. With the aid of this expression it is argued that in the case of (5-Pu, the SP is overtaken 
' by a variant of OP that does not break the time reversal symmetry as is confirmed by electronic structure calculations 
, within the LDA-I-U approach. In this context, there is a discussion of the nature of this exact expression for the 
OP and how it compares with earlier formulations 1^, 14 1. Finally, the corresponding OP multi-pole is discussed in 



, connection with existing d to f branching ratios experiments 13, Hi, where it appears through a sum rule. 
' In the most general version of LDA+U [l^l the HF correction enters with a Hartree (H) and exchange (X) term as 

^: 

Eh+ Ex^ [pacPbd - PadPbc) {ab\g\cd) , (1) 

abed 

where pab is one element of the density matrix for the i^^ shell, with dimension(4^ + 2) x {M + 2), (or 2[£] x 2[£] if we 
use the conventional notation, [£] = 2£ -I- 1) which acts as an occupation matrix. Here a is a combined label for the 
magnetic quantum number rua and the spin variable Sa- The interaction has the form 

{ab\g\cd) ^5{sa, Sc)5{sb, Sd)[ef 

X ^^'^c('^)(m„ mc)c('=)(mfc, m^). (2) 

kq 



where c^^^ are Gaunt coefficients and F^^^ are the Slater integrals of the screened Coulomb interaction, in this work we 
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stay with the convention and choose the latter on physical grounds. In the exchange term the spin Kronecker-deltas 
will allow for a non-diagonal spin interaction between the two density matrices p, giving rise to a spin-mixing. 

This method has been implemented in the full potential augmented plane wave (FP-APW) package Exciting 
17| . and the results above have been verified. A straight-forward density functional approach leads to an anti- 
ferromagnetic order with large spin and orbital moments, while switching on the LDA-I-U HF interaction, including 
spin-mixing terms, leads to a non-magnetic solution as displayed in Fig. (fTl) for a double counting of the type around 

□ 




U(eV) 

FIG. 1: Moments and exchange energies from LDA-I-U- AMF calculations of fee Pu (a=4.64A) within the FP-APW method. 
Basis set cutoffs corresponding to RmtGiuhx = 9, with muffin tin radius Rmt = 1.56A, and a Brillouin zone sampling of 864 
points for a two atom cell were used. The Slater parameters are summarized in two parameters U and J in the same way as 
in Ref. The spin (blue triangles) and orbital (red squares) moments are shown for a varying U but a constant J = 0.68 
eV, which corresponds to E^^^ = 53 meV, in the bottom part. Filled symbols indicate the LDA results. Also displayed are the 
5/ occupation number (black) as well as the w^^" (green) where the dashed green line indicates the corresponding saturation 
limit. At the top, the most significant terms of the exchange energy in Eq. (|6} are displayed as lines, with same color scheme 
as for the moments, with the tensor rank identifying each, given in the legend. 

To straight-forwardly analyze this non- magnetic solution is cumbersome due to the 4[£]'^ independent elements of 
the density matrix p. Instead we take the following approach. As a generalization of the fact that one gets the spin 
and orbital moment from the density matrix as e.g. (S) = Tr S p, we introduce the expectation value of double tensor 
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operators w'^p = Trv*^ p [l8, 3, 2^. In our basis the tensor operators take the form 



£ k e 
-■nib X m.a 



t^^(.,|t^K) = (-r-(^4^ ^ ;Jn;/. (3) 

Here we have use the so-caUed 3j-symbols (...) [2l[ and the same normahsation as in Ref. [l^. 

An attractive property of the multi-pole double tensors w'^p is their simple interpretations. As have been pointed 
out [3] the w'^^, with k even, are related through Wigner-Eckhart theorem to the A:-th multi-pole moment of the £ 
charge density while the w*''^, with even fc, are related to the multi-poles of the magnetization density. Finally the 
odd k correspond to currents, i.e. the tensors are related to the multi-pole moments of the current {p = 0) and the 
spin-current (p = 1). 

It is fruitful to view the introduction of w'^p as a transformation of the density matrix p to these double tensors. 
This transformation is one-to-one p^ . which is consistent with the fact that the number of parameters are kept. The 
sum of the rank of the tensor operators are X^fcl^] ^ SpH ~ which is equal to the number of independent 

(real and imaginary) components of p, since it is liermitian. So we keep the same information as in the density 
matrix, but now distributed on [k][p\ independent double tensors, which turns out to be very convenient. The inverse 
transformation is readily obtained by utilizing orthogonality properties of the 3j-symbols |2l| . 



kx 



X 

py 



-rric X ma 



Y.[pM-r-^( 4 ; X- (4) 



and 



As can be straight-forwardly shown, with the introduction of the double tensors w'^^' and the relation Eq. ([4]), both 
the direct and the exchange term can be put in simple forms, 

Eii = Y,F'-''h{e,k,k)w''° ■w''^ (5) 

k 

It 1 

k fci=0 p=0 

with 

n,„o^^(^SO" 

J(£, fc, fci) = /(£, fc, k{)^^^L^ 

where the {. . . } symbol is the 6j-symbol [2l|. This is a remarkably simple reformulation of the exchange energy. It is 
a straight-forward generalisation of the Stoner-likc exchange formulation of SP, with scalar products of quantities with 
themselves. The final expression is actually almost identical to an expression derived by Racah for a P' configuration 
[20| . but is more general in that it is valid for any non-integer occupation of the £-shell. The calculated interaction 
strengths from Eq. ([8|) become even simpler when going from Slater parameters to a certain linear combination of 
them called Racah parameters, as many values then vanish. The relation between the two types of parameters is given 
for the case of /-shell in Table |T] and the corresponding transformed interaction strengths are tabulated in Table |TT1 

The independent terms in this exchange energy expansion have simple physical meanings, for instance the term 
with w'^^ is the Stoner-like SP, while the terms involving w*^"^^ with k\ — 2,4,6 are responsible for an intra-atomic 
non-coUinear 

spin polarisation B. In particular , we have an explicit expression for orbital polarization, or Hund's 
second rule, for a general non-integer system. From Eq. ([S]) and with the observation that orbital moments are 
(L) = £w^°, we get a generalized OP term for / states as 
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TABLE I: The expansion coefficients c{k,k') of tlie Racah parameters for the /-shell, = c(k, k')F^'^ \ are given. We can 
note the relations with the standard U and J parameters: -B*-"' = U — J and E^^'' = | J. 
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TABLE IL The interaction strenghts J(3,fe, fei) from the multi-pole expansion of the exchange energy in terms of Racah 
parameters, E'^''\ = - X;^^^^ -E'*' J(£, fc, fci)w*^if ■ w*^!" . 
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This expression is very close in resemblance to the expression suggested by Brooks 13, [lj| (OP-B), which in our 
terminology looks like 



E^-^^-^yiSwl^,)- . (10) 



p 

We note that there are three corrections to this simplified OP-B formulation. Firstly, there is a contribution from 
£^(0) = U — J too. Secondly, for /-states the exact contribution proportional to E'^^^ is 33/112 « 0.295 instead of 
1/4. These two corrections can in principle be compensated by using a somewhat larger effective i?^"^-* parameter in 
OP-B. However thirdly, in OP-B only the z components are included (tensor components 00). This is due to that spin 
mixing are neglected. This will be found to be a more severe restriction, since in general we can get an independent 
contribution from each of the three spin components. All these three corrections lead to an effectively stronger OP 
than what was originally suggested by Brooks. In addition wc observe that the OP is essentially two different terms; 
one that favours broken time reversal symmetry states with w^*^ ^ (OP-odd), while the second term involving w^^ 
does not break time reversal symmetry (OP-cven). 

There have been other attempts to improve on the OP-B formula (23 . [25j by starting from an integer occupation in 
the same spirit as the original work [T^]. However, since there are uncertainties in this limit what should be attributed 
to OP, the resulting formulations usually involve other terms of Eq. [6] and there have never been a general description 
of OP for a / shell. The only approach we are aware of for non-integer systems also failed to get an expression for / 
systems 
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Let us return to the LDA+U calculations of Fig. ([T]) where the different energy contributions arising from the 
different tensor products are given. Here the two angular momenta kp of the double tensor are coupled into a third 



as discussed in Rcfs. 18|, |20|, giving rise to an irreducible tensor w'^'^''. As can be seen only a few tensors have 
any significant contribution to the total exchange energy. They arc w°"° (total 5/ occupation), w"^^ (SP) and w^^'^ 
(OP-even). The OP-odd term w^"^ is almost detectable for low U values. It is worth noting that the total exchange 
energy calculated by Eqs. ([T|) and ^ are indistinguishable. It is evident from the graph that the OP-cven term 
takes over the SP exchange energy when increasing the effective Coulomb interaction U. In fact, as the AMF double 
counting corresponds to neglecting the contributions to the HF exchange from w™" and w°^^, we observe that the 
OP-even term solely determines the HF exchange for all values of Ul By studying the calculated observables as 
displayed in bottom part of Fig. Ilj it is worth stressing that although there is a small increase in the 5/ occupancy 
from 5.2-5.5, it has little influence, in contrast to what has earlier been assumed 0,01 • It is the steep increase in the 
magnitude of w^^'^ that stabilises the non-magnetic state. 
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The identified broken symmetry of the calculated state is quite intruiging. The order parameter w^^° corresponds 
to that the three components of the spin currents orbit around their different spin quantization axes with equal 
magnitudes. This leads to a time reversal invariant, scalar, order parameter, and since it arises from spin currents it 
is a quantity difficult to observe directly in experiments. 

Recently there have, however, been reports on measurements on the branching ratio for the d to 5 / tr ansition 
for several actinide systems [13, [lH, from which values of w^^*' can be obtained through a sum rule [2^. These 
measurements report very large values, not least for a-Pu. In the subsequent discussion they attribute this to the 
strong spin-orbit coupling which brings the 5/ states close to a jj coupling scheme. In the light of our finding we 
would like to alter that analysis slightly. While the spin-orbit coupling is important in the actinides, it is not strong 
enough to bring the 5/ states into a jj-limit by itself. In fact without the HF term of Eq. ([1]), i.e. in the LDA limit, 
we calculate a spin-orbit-only value of —2.4 while in the presence of the HF term wc get enhanced values varying 
between —4.4 and —7.2, as seen in Fig. The values for large U parameter are close to saturation, as indicated 
in Fig. (Ill), which would correspond to — Iw^"". These values should be compared to the measured value of —5.1 
for Pu in its a-phase, assuming a configuration [T^]. We notice that the exchange term is essential to bring the 
calculated w^^*' to the same magnitude as the experimental value. This leads to the conclusion that there is a strong 
competition between different exchange channels in the actinides, where the spin-orbit coupling plays a role since it 
favours the OP-even channel over the SP channel. This in accordance with calculations on other actinide systems, 
where we have found that the w^^° always have a large contribution, even for magnetic systems [lit . 

As a summary, we conclude that while all calculations essentially involve the dominant SP and OP contributions 
to exchange, the main difference between the calculations leading to a magnetic state [3, [^ and the ones leading to 
a non-magnetic state [3, [3| is their treatments of the OP term, the former utilizes OP-B without any spin-mixing, 
while the latter uses the correct OP-even including the spin-mixing contribution. This symmetry broken state, with 
non-zero spin currents, has a surprisingly simple structure since all non-trivial exchange energy goes into the OP-cven 
channel. Since the exchange energy usually is larger than the correlation energy it is likely that this state has many 
similarities with the ground state of ^-Pu. 

The derived multi-pole expansion of the exchange energy is general and valid for all types of open shells. A more 
general study of the relevant exchange terms for other / and d systems will be published elsewhere [3 • 

This work illustrates a large advantage with the multi-pole expansion of the exchange energy of Eq. ([6]) - it brings 
forward the physically important exchange channels in a simple way. It is clear that most part of the terms in the 
expansion has little or no contribution. Hence, in the future it is of great interest to perform similar multi-pole 
expansions on the terms relevant for studies of correlation within the / shell, e.g. the Green's function and the self- 
energy. In particular, in order to better understand S-Pu and its intriguing electronic structure, one ought to study 
the effect of correlations on the spin currents of the OP-cven state, by means of e.g. DMFT calculations together with 
a similar analysis as performed here. 
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